Driving-induced crossover: from classical criticality to self-organized criticality 
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We propose a spin model with quenched disorder which exhibits in slow driving two drastically 
different types of critical nonequilibrium steady states. One of them corresponds to classical critical- 
ity requiring fine-tuning of the disorder. The other is a self-organized criticality which is insensitive 
to disorder. The crossover between the two types of criticality is determined by the mode of driving. 
As one moves from "soft" to "hard" driving the universality class of the critical point changes from 
a classical order-disorder to a quenched Edwards- Wilkinson universality class. The model is viewed 
as prototypical for a broad class of physical phenomena ranging from magnetism to earthquakes. 
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The study of criticality in externally driven inhomoge- 
neous systems has attracted much attention in the last 
two decades [Ij. Sueh systems exhibit rate independent 
dissipation and are widely used to model hysteretic phe- 
nomena and intermittency associated with magnetism, 
superconductivity, porous flow, fracture, friction, plas- 
ticity and martensitic phase transitions [3]. It has been 
realized that in most cases scaling emerges as an inter- 
play between quenched disorder, extremal dynamics, and 
quasi-static driving [2J|. Within this general framework 
the theoretical work has been mostly focussed on two 
types of models. In models of the Random Field Ising 
(RFIM) type, the critical behavior requires fine tuning 
of the amount of disorder r and the intermittent events 
(avalanches) are scale- free only at a certain r — ro [l|. 
An alternative approach links criticality to a pinning- 
dcpinning (PD) transition where the disorder r is an ir- 
relevant parameter [3[ . It has been established that criti- 
cality in the first class of models is classical, as in second 
order phase transitions [l|, while in the second class it is 
self tuning in the sense that infinitely slow driving brings 
the system automatically on a critical manifold [J, la, ly] . 

The two approaches are fundamentally different. The 
first model describes regimes with dominating nucleation, 
while the second one deals exclusively with propagation. 
It is then not surprising that the resulting criticality is 
different. In the RFIM the emerging scaling has been 
explained by the existence of a classical critical point 
of the order-disorder (OD) type. On the contrary, in 
PD theory one encounters a range of universality classes 
none of which can be formally reduced to OD. A rela- 
tion between the two types of critical phenomena has 
been, however, anticipated. Previous work has shown 
that the presence of a nonlocal demagnetizing field of 
antiferromagnetic nature (as in soft magnets) can self 
tune the RFIM to display front-propagation critical ex- 
ponents [71, Isl]. In the PD framework similar 'self-tuning', 
often interpreted as self organized criticality (SOC), is 



obtained if the system is driven through a 'weak spring' 
which provides an explicit feedback mechanism [9|. In 
this letter we use these insights to develop the first uni- 
fying model which displays a crossover between the OD 
and quenched Edwards- Wilkinson (QEW) universality 
classes. We show that such crossover can be achieved 
experimentally by modifying the properties of the exter- 
nal driving. Since the QEW model is equivalent to the 
Oslo rice pile model [1^ and is therefore paradigmatic for 
SOC, we are essentially dealing here with a fundamental 
relation between OD and SOC. 

We base our model on the observation that solids can 
be deformed either by applying a force (soft device), o r by 
controlling surface displacements (hard device) [ISJ, [iJ] . 
Both driving mechanisms, soft and hard, can be handled 
simultaneously if the system has a finite 'elasticity'. To 
introduce this effect in a spin setting, consider a proto- 
typical model which we call RSSM (Random Soft-Spin 
jlSj or Snap-Spring [l6|, ll7| Model). The main differ- 
ence between the RSSM and its predecessor RFIM is 
the finite curvature of the energy wells and the pres- 
ence of elastic barriers. While the role of the softness 
of the spins is known to be secondary in equilibrium and 
in the case of soft driving, where it can be accounted 
for through an appropriate 'dressing' of the underlying 
hard-spin model [l5| , the curvature of the wells becomes 
crucial in the case with hard driving where it plays the 
role of a regularization of the otherwise degenerate prob- 
lem 1^ . Consider a set of N bi-stable units located on 
the nodes i = l,2,...7Vofa cubic lattice with linear size 
L = N^/^. The state of each snap-spring is characterized 
by a continuous scalar order parameter e^ measuring the 
local 'strain'. The energy of the system is 
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where fi{ei) is a double- well potential and K = {Kij} 
is the interaction kernel with sufficient rate of decay to 



ensure convergence in the thermodynamic Umit. In each 
well (one defined for e^ < and the other for e^ > 0) we 
use the approximation fi{ei) = ^{ci — s,;)^ — h'^Ci, where 
Si = ±1 is a spin variable and {/i^} are random numbers 
representing quenched disorder. The system is loaded 
through an 'elastic' device with the energy cpd = |(e — e)^, 
where e = (1/-/V) J2i ^i is the average strain of the sys- 
tem of snap-springs, c is the stiffness of the loading 
device and e is the control parameter (global 'strain'). 
One obtains a hard device in the limit c —>■ oo, and a 
soft device in the limit c — > 0, e — > oo with the stress 
a ~ ce fixed; the system will be driven quasi-statically 
by changing the control parameter e, if c is finite, and 
cr, if c = 0. We neglect thermal fluctuations (T = 0) 
and assume that the harmonic variables e^ relax 'instan- 
taneously' and can be adiabatically eliminated. By mini- 
mizing the total energy (t>+(j)d with respect to e^ we obtain 
Ci = e + J2f=i {J'lj - k/N) Sj + h, where J = (1 -f^ K)-i 
is the effective interaction between the spin variables Si. 
We impose periodic boundary conditions, meaning that 
koo = X)i Jij does not depend on i, and use the notations 
e = e[c/coo(cfcoo -f 1)"^], k = koo[ckryo{^ + ck^)-\ and 

f^i — 2.^ j ij j ■ 

Because of the absence of thermal fluctuations and to 
the separation of time scales between overdamped relax- 
ation and driving, the system remains on a metastable 
branch {ei({sfc};e)} corresponding to a particular local 
minimum of the total energy until the latter ceases to 
be stable [13|- When the instability condition SiCi < 
is reached for some i the system jumps (through an 
avalanche) to another locally stable branch character- 
ized by a different spin configuration {s[}. We increase 
e by driving the system from an initial stable configura- 
tion with {si = —1} and assume that avalanches prop- 
agate (at constant e) with synchronous dynamics. For 
simplicity, we consider only nearest-neighbor interactions 
with Jii = Jo, Jij ~ Ji, and we set Jo = ^i == 1 so 
that koo — 7. The renormalized disorder variables {hi\ 
are drawn from a Gaussian distribution with zero mean 
and standard deviation r. Under these assumptions, the 
RSSM becomes formally equivalent to a nonlocal aug- 
mentation of the classical RFIM with demagnetizing fac- 
tor A: [i|,[i,0,[i,[il. 

At A: = we expectedly observe an OD transition 
at To — 2.2 which separates a 'popping' (POP) regime 
(r > To) where all the avalanches are small from a 'snap- 
ping' (SNAP) regime (r < To) where an infinite avalanche 
sweeps most of the system [l|, |l9[ . The stress-strain curve 
cr(e) = k^{e — kN^^ J2i ^i) is continuous in the POP 
regime and displays a macroscopic discontinuity of the 
strain in the SNAP regime [see Fig. [1] . The strain dis- 
continuity associated with the nucleation of the infinite 
avalanche occurs at a 'nucleation' stress (j„ which de- 
creases with r as indicated in Fig. [51 

For non-zero fc, the POP regime remains essentially 
unaltered. In contrast, stiffness has a non-trivial effect 




FIG. 1: Stress-strain curves for A; = (dotted line), k = 
0.5 (dashed line), and k = koo (continuous line) for r — 1.5 
in a system with L = 64; cTp and (t„ are propagation and 
nucleation thresholds, respectively. 
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FIG. 2: Dependence on disorder of the averaged nucleation 
and propagation thresholds: {a„) (squares), {ap) a.t k — koa 
(circles), and {upd) for an interface [100] (triangles). Insets: 
cross sections of the 3D system showing typical transforma- 
tion domains for k — koo- Darker and lighter colors indi- 
cate transformed regions (s — +1) and untransformed regions 
(s = -l). 



over the SNAP regime observed at low disorders. In 
this case, a compact domain reminiscent of the infinite 
avalanche (as in the k = case) nucleates at the stress 
(T„. The stress relaxes during the avalanche growth to 
satisfy the global driving constraint. When k is below 
a certain threshold kp{r), the SNAP behavior remains 
qualitatively as in the fc = case because the nonlocal 
constraint is still soft. In contrast, when k > kp(r), the 
first stress drop prevents the spanning avalanche from 
growing [Fig. [T] . The transformation induced by the sub- 
sequent increase of e proceeds through the intermittent 
growth of the previously nucleated domain with untrans- 
formed system acting as a disordered background. The 
propagation is accompanied by stress oscillations around 
the propagation threshold CTp which remains stable dur- 
ing the whole yielding process. The average of Up over 
disorder coincides with the depinning stress apd for a 
flat [100] interface artificially introduced and driven as 
in Refs. [20|, l2]J. This is a clear evidence that systems 
with k > kp{r) reach a front-propagation regime which 
self-tunes exactly around the PD point. The surface mor- 
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FIG. 3: Scaling collapse of Nsik,r,L) (a,b) and D{S;k,r,L) 
(c,d) according to the hypotheses in Table IH Different sym- 
bols correspond to different system sizes as indicated in the 
legend. Data in plots to the left (right) correspond to r = 1.5 
(r = To). Data in (c) correspond to wL''" = 66 ; in (d) to 
kL''° = 400. The lower scaling collapse in (d) is obtained with 
exponent r, while the upper one corresponds to To; dashed 
lines indicate the power-laws expected in the thermodynamic 
limit. 



phology of the growing domain is faceted at low disorders 
(r < 0.8 for systems with L = 100) and is sclf-affinc 
at intermediate disorders. According to previous studies 
221 . the faceted morphology is unstable giving rise 



[211, 

to self-afSne morphology in the thermodynamic limit for 
any finite r. Given the self-affine character of the bound- 
ary of the nucleated domain in the case of intermediate 
disorders and the short-range character of the involved 
interactions, one expects to observe critical scaling of the 
QEW (quenched Edwards- Wilkinson) universality class 

We check the validity of this assumption by analyz- 
ing the scaling properties of the number of spanning 
avalanches Ns{k,r,L) and the distribution D{S;k,r,L) 
of avalanche sizes S. Following the arguments given in 
Rcfs. [19[, we consider contributions to Ns of avalanches 
spanning along 1 or 2 dimensions only. Table [J sum- 
marizes the action of the Renormalization Group (RG) 
transformation on the parameters of the system and the 
resulting scaling hypotheses. The scaling variables mea- 
suring the distance to the OD critical point are u = r — Vo 
and k. The distance to the QEW critical manifold is 
w ~ k — kpir). 

Scaling collapse for Ns{k,r,L) close to the line kpir), 
shown in Fig. [H^a), confirms the validity of our scaling 
hypothesis. We have checked that QEW scaling persists 
over a finite interval of disorders r < Tq. QEW scal- 
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FIG. 4: (a) Schematic phase diagram for the RSSM (thermo- 
dynamic limit). Criticality of the OD type is associated with 
the point r = ro,k = Q. Criticality of the QEW type is ex- 
pected on the line kp(r) separating SNAP and POP regimes; 
symbols display estimations of fep(r) from scaling collapses. 
(b) Schematic RG flow for the RSSM model. Separatrix 1 
going from the neighborhood of the OD flxed point towards 
the QEW flxed point indicates the QEW critical manifold. 
The RG-flow towards SNAP and POP regimes is indicated 
by arrows 2 and 3, respectively. Lines 4 and 4' correspond to 
systems which display SOC as an intermediate asymptotics 
(QEW exponents with supercritical or subcritical cutoffs for 
k < fcp(r) and k > kp{r), respectively). 



ing fails for disorders close to Tq (above r ^ 2 for our 
system sizes) due to crossover to the OD critical regime. 
Fig. ^h) displays the scaling collapse at u = which gen- 
erates the exponents po and 9o given in TablelD The value 
of 6o which, in contrast to 9q is positive, agrees with pre- 
vious estimates for the RFIM [19|. Fig.jSjc) presents the 
scaling collapse of D{S; fc, r, L) for r = 1.5 corresponding 
to a particular section of the function D{SL^'^'','wL'^'>) 
at constant wL''« . The resulting exponents dq and Tq 
listed in Table |T] do not depend on the selected section 
and are in agreement with previous estimates from front- 
propagation models \4, lD, l21| . This confirms once again 
that the propagation regime is of the QEW class; the 
analysis of statistics of durations, omitted here, also sup- 
ports this interpretation. The value of dq, giving the di- 
mension of the avalanches in the propagation regime, is 
consistent with d—1, which confirms the self-affine mor- 
phology of the propagating domain boundary [2l| . The 
distribution D{S] k, r, L) at r = tq scales with exponents 
previously reported for the OD universality class when k 
is very small. The exponent Tq displays a considerable 
crossover to QEW when the stiffness becomes relatively 
large. For instance, Fig. Eljd) shows that the scaling of 
D{S; k, r, L) at kLP° = 400 is better with exponent Tq 
than with Tq. 

In Fig. mja) we present the phase diagram in the (r, k) 
plane showing stable SNAP and POP regimes separated 
by the QEW fine u; = 0. This fine describes QEW be- 
havior and it ends in a point corresponding to the OD 
regime. The variety of observed non-equilibrium steady 
states can be explained if one assumes the existence of 
four fixed points for the RG flow, schematically depicted 



TABLE I: Scaling of the relevant quantities under a RG transformation with blocking parameter b and values of the associated 
critical exponents (close to OD (subindex 'o') and to QEW (subindex 'g')). 



OD 



QEW 



RG 



exponents 



RG 



exponents 



System size, L L{b) — b L 
k{b) = 6"° A: 



Stiffness, k 
Disorder, r 
Size, S 

Ns{k,r,L) 
D{S;k,r,L) 



Ns{k,r,L) 
D{S;k,r,L) 



S{b) ^b-"^"S 
Ns{b) = b-'>"Ns 
D{b) =6"°''°£) 



Po 

do 



1.3 ±0.3 



1.2 ±0.1 (fW]) 
2.78 ± 0.05 ([19|) 
60 = 0.10 ±0.02 ([19|) 
To = 1.6 ±0.06 (UJ) 



L(6) = b-'L 
w{b) = b'''"w 
Irrelevant 
Sib) = b"^''S 
Ns{b) ^ b-"" N, 
D{b) = V'^-'D 



Pq =0.8 ±0.2 

dg = 2.0 ±0.1 
Oq = -0.2 ± 0.06 
r„ = 1.3 ±0.06 



L"''NsikLP'',ku- 



') 



ku 



L"" NsiwL"") 



in Fig. mjb) . The OD regime is associated with a fuUy re- 
pulsive critical point which can be reached only by tuning 
all four parameters: a = Oc, r = ro, k = 0, and L~^ = 0. 
In contrast, QEW is a saddle point with a stable mani- 
fold which governs the large scale behavior of the systems 
with r < To, CT = CFpdir), k = kp{r), and L~^ = 0. As we 
have seen the condition a = apd{r) is reached automat- 
ically during the self-organized propagation regime with 
k > kp(r)\ the corresponding systems lay on the criti- 
cal manifold connecting OD and QEW points. The large 
scale behavior for systems located away from the criti- 
cal manifold is governed either by SNAP (arrow 2) or by 
POP (arrow 3) fixed points which are trivial attractors 
in the RG sense. While the self-organized propagation 
regime is strictly critical only for k = kp(r), our numeri- 
cal simulations show that the system exhibits truncated 
power law scaling with QEW exponents in a broad range 
of parameters around the line kp(r) (arrows 4 and 4' in 
Fig. m^b)). This observation suggests that SOC can be 
viewed as an intermediate asymptotics. In general, since 
in real systems stiffness is finite and disorder is generic, 
we anticipate the power law structure of critical fluctu- 
ations to be in most cases of the PD rather than of the 

OD type. 
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